arXiv:gr-qc/0211056v3 24 Dec 2002 


EPHOU 02-007 


Analytic Evaluation of the Decay Rate for an Accelerated Proton 
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We evaluate the decay rate of the uniformly accelerated proton. We obtain an analytic expression 
for inverse beta decay process caused by the acceleration. We evaluate the decay rate both from 
the inertial frame and from the accelerated frame where we should consider thermal radiation by 
Unruh effect. We explicitly check that the decay rates obtained in both frame coincide with each 
other. 
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I. INTRODUCTION 

In 1974 Hawking announced a celebrated result that 
the quantum thermal radiance occurs from black hole [Q . 
which is the most interesting result of the quantum field 
theory in curved space-time. Since then, the study of 
the field theory has been of much interest. It appeared 
that the quantum radiation is not special phenomenon 
in strong gravitation but analogous effect exists even in 
the flat space-time [|). The effect is caused by the fact 
that there is no vacuum which is invariant under the 
general coordinate transformations although the vacuum 
is invariant under the special coordinate transformations. 
A typical example is the accelerated observer which feels 
thermal radiation of the particles depending on the rate 
of acceleration. The effects in the accelerated system are 
not hard to be analyzed and is more tractable than the 
purely gravitational effects. 

This phenomenon in flat space-time was shown by 
Fulling-Davies-Unruh in 1973, which is now called Un¬ 
ruh effect H [|. A vacuum in rest frame is not a vac¬ 
uum in accelerated systems. Actually, the vacuum in 
rest frames is a thermal bath in the accelerated frames 
Q. This means that the definition of particle in quan¬ 
tum field theory depends on the observer and a strongly 
curved space-time by the acceleration corresponds to a 
strong gravitation. 

Until now, many phenomena which the acceleration 
gives the interesting influence have been shown in the 
cosmological setting |j. There is only one experiment 
about Unruh effect by Chen and Tajima |7j. They found 
that an electron which is accelerated by electronic field 
quivered. They showed that the random absorption of 
quanta from the FDU thermal bath cause this quivered. 
However, the experiment of decay is not existing. 

More recently, the decay of accelerated particles was 
studied. In 1997 R. Muller [H estimated the decay rates 
of the accelerated particle [||. However, it was shown 
that decay rate is very small in our experimental setup. 
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who realized that the decay rate can be used as a 
”theoretical check ” of Unruh effect. The most interest¬ 
ing decay is in the case where we can deal with a particle 
which doesn’t decay in inertial frame but can in accel¬ 
erated frame. A typical example is for the proton. In 
inertial frame the accelerated proton can decay, and we 
can write the process as 

r res t; P + —► n + e + + v e . 

In the accelerated frame, the observer takes the same 
acceleration as proton then the proton is stable but his 
space is in a thermal bath of particles. So the proton 
should absorb e~ and v, and the processes become the 
inverse [3 decays as 

F a c c ; P + + e~ -> n + v, 
p + + V —> n + e + , 
p + + e~ + V —> n. 

Both of the cross sections should agree with each other 
because they are the probabilities of the identical physi¬ 
cal phenomenon and it is independ on the frames we use. 
However, the calculation is very complicating because the 
solutions of the Dirac equation in accelerated frame be¬ 
come a special function even if we represent protons by 
the classical current and restrict to tree level. There¬ 
fore, they restricted the system to two-dimensions. They 
could represent the decay rate analytically for Minkowski 
frame in the limit where neutrino is massless. However, 
the decay rate for the accelerated frame could not be 
evaluated analytically. Therefore, they showed the agree¬ 
ment numerically for massless neutrino. This calculation 
shows that the Unruh effect is inevitable for the accel¬ 
erated frame. An interesting point of their approach is 
that they use a decay rate, true physical observables, for 
showing the existence of Unruh effect although the effect 
itself was established earlier (see Ref. |J for Review). 

They used two dimensional setting and treated neu¬ 
trino for massless field for simplicity. However, it is de¬ 
sirable if we can show that decay rates calculated in both 
frame are identical even for massive neutrino in four di¬ 
mensions because massiveness of neutrino is now proved 
exp er iment ally. 

The main aim of this paper is to show that we can com¬ 
plete their calculation by showing that the decay rate of 
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proton is independent of the frame not by numerical com¬ 
putation but by analytical computation. Moreover, We 
will perform the evaluation in four dimensional setting 
and treat neutrino a massive field. It is interesting that 
the decay rate can be analyzed analytically. 

In Section II we will calculate the cross section of (3 
decay in the inertial frame. We will show that the decay 
rate can be obtained analytically in terms of a function 
which is an analog of Meijer’s G-function of two variables. 

In Section III, we will perform the calculation in accel¬ 
erating frame. We will be able to check that the resulting 
function is identical to the one obtained in inertial frame. 

Section IV is devoted to the discussions. 

In Appendix B, we list the explicit form of the function 
appeared in our main result. 

II. INERTIAL FRAME 


where q is a small coupling constant and u M is the four- 
velocity u^ = (a, 0,0,0) in Rindler coordinates. To deal 
with the proton-decay, nucleons | n) and protons | p) are 
represented as excited and unexcited states of the nu¬ 
cleon, respectively. For Hamiltonian H , neutron and pro¬ 
ton mass m n and m p , respectively, we set 

H\n) = m n \n), H\p) = m p \p). (4) 

We replace q in last current by the Hermitian monopole 
q(r) = e^me- i6 \ ( 5 ) 

where r = v/a is proper time of the proton. 

Then the current four-vector is 

j M = q(r)u tl S(x 1 )5(x 2 )5(u — a -1 ). (6) 

This is clear that the particle is on the hyperbola. 


In this section, we analysis the (3 decay of the accel¬ 
erated proton in the inertial frame. In this frame, the 
accelerated proton decay resulting from the acceleration. 

In the flat space-time, the solutions of the Dirac equa¬ 
tion are simple but the calculation becomes complicate 
because of the vector current is on the hyperbola. We 
will mainly follow the notation appeared in Ref. 00 


A. Accelerated Proton current 


There may be several methods to represent the accel¬ 
eration. One way is the followings. We represent the 
proton as a classical current. The position of the current 
should be on the world line of the accelerated particle, 
“hyperbola”. To do this, we introduce the Rindler co¬ 
ordinates. Rindler coordinates correspond to the world 
line of the uniformly accelerated observers. The inertial 
coordinates are shown by (a: 0 , cc 1 , x 2 , x 3 ) and the Rindler 
coordinates which show acceleration along the 3-axis are 
(v, x 1 , x 2 , u) with 0 < u < oo and — oo < v < oo. Two 
coordinates are related by 


x° 

= ztsinhu, 

1 


X 

= x , 

2 

2 

X 

= x , 

x 3 

= u cosh v 


(1) 


Then the line element is described in the Rindler coordi¬ 
nates by 


ds 2 = u 2 dv 2 - (dx 1 ) 2 - (dx 2 ) 2 - du 2 . (2) 


If some particle uniformly accelerated with a proper ac¬ 
celeration a, then u = a -1 = const, is its world line. As a 
result, the protons which are accelerated to the axis of x 3 
on x 1 = 0, x 2 = 0 are represented as following classical 
current 


( 3 ) 


B. Fermionic field quantization 

For electrons and neutrinos, we write the fermionic 
fields satisfying the Dirac equation (*7 M <9 M — = 0 

as 


41 (x) = 


CT —± ' 


d 3 k 


b ka ip£ u) (x) + dl a tpi_ u J (x) 


( 7 ) 

where x = (a; 0 ,x x ,x 2 ,x 3 ), k = (w, k 1 , k 2 , fc 3 ), x = 
(cc 1 , x 2 , x 3 ), k = (k 1 , k 2 , k 3 ). bka and d^ kcr are annihila¬ 
tion and creation operators of fermions and antifermions, 
respectively, k and a are momentum and polarization, 
respectively, and V’k ^ are positive and negative 

frequency solutions of the Dirac equation. 

By solving the Dirac equation we obtain the orthonor¬ 
mal mode solutions [p^|, 


v-£“ 


'(*) = 




(2tt)' 


,(±V 


(k), 


where 


<4 ±w) (k) 


ku'y^ ± m 

=u, 

y/u)(u) ± m) 


and 



-1 ■ 


’O' 


0 


1 

U+ = 

0 

, U- = 

0 


.0. 


.0. 


( 8 ) 

( 9 ) 


( 10 ) 


We choose the traditional inner product of the form 


JT. 


= <5 3 (k - k')5 aa i6± UJ ± UJ ' (11) 


= qu fJ 'S(x 1 )S(x 2 )S(u — a : ), 
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and the solutions are normalized by this definition, where 
ip = ip' 7 0 , dT, p = n p dT,, n p is a unit vector orthogonal 
to £ and we have chosen £ to be the hypersurface of con¬ 
stant x°. The creation and annihilation operators satisfy 
the equal-time anti-commutation relations: 

{ka,bt r<7l } = {d klJ .dl a ,} =d 3 (k-k')S IJI7 , 1 

^ ^kcr , bk'tj' j" — d]tcr , dk'a' J" 

= = {&k<x,dE'} = 0. (12) 

Using electron tfe and neutrino d® fields, we write the 
Fermi action by the nucleon current (2.6) 


Si= / d^xyf—g)^ v + 


The way of deriving the cross section is normal. 

First, the vacuum transition amplitude of the proton 
decay is written by 


A p ^ n = (n\ ® (ef ,^ ( 7 JS' / |0) <g> | p). (14) 


It is straightforward to compute A p ~* n for the Fermi ac¬ 
tion Si and we obtain 


/ OO 

d{x°)d(x 3 ) 

-OO 


JAmr, 


sja 2 (x 0 ) 2 + 1 


(13) 


xS(x 3 - \/(x 0 ) 2 + a~ 2 )(e]t a , r'^ CT j4' l/ 7' 1 'I' e |0),(15) 


where the first term is used by inverse f3 decay. 

You can see these general formation on any book of 
Quantum Field Theory in the inertial frame. We are 
ready to start to calculate the cross section of the accel¬ 
erated proton in the inertial frame. 


where A?ti = m n — m p , r = a -1 sinh _i (ax°) is the nu¬ 
cleon’s proper time, Gf = (p|<)(0)|n) is the Fermi con¬ 
stant. We can substitute the field and integrate by x 3 . 
The differential transition rate is 


C. Calculation of Cross Section 


d 6 V p ^ n 

d 3 k e d 3 k l 


= £ 5 > 


p->n\ 


(16) 


We are now going to calculate the cross section of the 
/3 decay using the current and field in the subsection A. 


We obtain it by integration of proper times t± and r 2 


J 


dPvf: 


g 2 f 


dT\dT2 lin'd 


:££[ 5 £"- 




— ( A-UJu ) v ^ 

C ’l U-a e 


d 3 k e d 3 k„ (27t) 6 J_ c 

x expi [Ato(ti — t 2 ) + a~ 1 (u> v + w e )(sinhaTi — sinhar 2 ) — a~ l (k 3 + fcg)(coshari — cosh ar 2 )] .(17) 


Now in order to calculate the spin sums, we use the fol¬ 
lowing standard formula: 


EE 


u^r 2U ^ 


a a cr p 


= Tr 


(±^),-(±o;)^ 0 r t 7 0 ^ U { ± u) vf> ±u) 


r iE u fe"'^r'7' 


&/3 p 




(18) 


where a and (3 represent e or v. 

The completeness relations can be written as 

E u ^ W “ ) ( k «)^ W “ ) ( k «) = (19) 


and we introduce s and £ by 


Tl = S+ 2 ’ T 2 =S “ 2 - 


( 20 ) 


By using the spin sum 


EE [ 4 EVEE 


— (A-UJv) V (-^e) 


<Te Ov 

1 


LU e U>, 


■ (w e w„ + k 3 k 3 ) cosh 2 as 


— (0 o e k 3 + k 3 LU,y) sinh 2 as — mem^ cosh a£] 
+(odcl function of k * and fc^,) 


( 21 ) 


and the change of variables 


k' a = k^, k' 2 a = kl,, k ,3 a = — u a sinh as + k 3 a cosh as, 

( 22 ) 
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we can obtain the differential transition rate 


g 2 f 


1 f°° <*>' +oj' 

_ / ^ e »[AK+ “ e 2a "“ sinh^] 

J — oo 


i d 6 pr* n 

T d 3 k' e d 3 k' v (2tt) 6 u' e 

x (w'wj, + k! e k' v — m e m v cosha£), 

(23) 

where T = ds. To integrate by £ we redefine new 
variable by 

A = e^. (24) 

And we use the following notations 

r k 'a ~ ~ rn a ~ A to 

fc a = —, ui a = —, m a =-, A m= -. (25) 

a a a a 

Then the cross section is given in the form: 


The integral of A can be readily expressed as modified 
Bessel to find 


o2„5/-o2 r oo 

nr n = ~ ~ / d 3 k e d 3 K 


n e e nAm J 0 
1 ffleffl,. 


K, 


2iAm 


K. 


2iAm-\-2 


+K, 


2iAm—2 


^ 2 (ljJ e + Wi,)j 

^ 2 (w e + w,,)j 

^ 2 (w e + Wj/)^ | . 


(27) 


However, this form is difficult to integrate by k e and 
Therefore, we use the integration formula of modified 
Bessel 


rJT" = A / -or* 

r,5Q2 r°° 


d 3 k e d 3 k„ 


2 5 7T 6 


dA e *( S e+W„)(A-^) 


> — oo 


xA 


2iAm —1 


c^eo;^ - -m e m u 


+ ^2 ) 


(26) 


1 /* r] q / Z \ 2s+/^ 

■ (28) 


It is not hard to show this formula in complex s plane by 
picking the residues of T(— s) and T(— s—p) (see Appendix 
A). To evaluate Eq. (E7h we use this formals and after 
some simple shift of variable we have 


n 5/~i2 
pp—m 'gt p 


7 r 6 e 7rAm Jo 


d 3 k e d 3 k v [ ( \/k 2 + rh 2 + \/ k 2 + m 2 


>c i 


2s-\-2iAm 


_—-_ 1 jji rn f - —-_ -—-_ "j 

r(—s)r (—s — 2iAm) + - J'A \ r(—5 + 1)T(—5 — 2iAm — 1) + r(—5 — l)r (—s — 2iAm + 1) \ 

2 co e uj v l J 


,(29) 


where the contour C\ must be chose so that all poles of T(—s)T(— s — p) are picked up and must be selected as k 
integration don’t infinity. But, this form of integration is still not simple for k e and 
In order to perform the integration, we use the following expansion formula 


(a i p\z __ f dt T( t)T(t z ) . t+Z p t 

{A + B) ~ J a Si r(- 2 ) 4 B 


(30) 


to integrate easily by k e and k v , where the contour C 2 is the path separating the poles of T(—t) from those of T(t — z) 
(see Appendix A). 

Now we can rewrite Eq. (129) as 


pp^n _ 


a 5 G 2 F 


ds 


dt 


(ml) s (miy 


2 ^7r 2 Jc s Jet r(— s — t + 3)F(— s — t + ^) 


r(-s -1 + iAffl + 3) r(-s)r(-*)r(-s + 2 )T(-t + 2 ) 


+Re |r(-s — t + iAm + 2)T(—s — t — iAm + 4)| T(—s + ^)T(—t + ^)T(—s + ^)T(-t + ^) 


(31) 
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where the contour C s and Ct is the path which picks up 
all poles of Gamma functions in s and t complex planes, 
respectively. 

This is the two-dimensional analog of Meijer’s G- 
function | ji~3| . The explicit form can be obtained by eval¬ 
uating contour integral. We list the results in the Ap¬ 
pendix B. 


then we find that the Dirac equation takes the form 

1 


6 lJ udi(udjX 1 ) = 


m 2 u 2 + - — Lo 
4 


5 lJ udi (udjX 2 ) = m 2 u 2 T - — w 


2 x'i - «<Dcr 3 X2,(40a) 
2 X 2 - iwcr 3 xi,( 4 0 b) 


III. ACCELERATED FRAME 

In this section, we are going to analyze the same phys¬ 
ical phenomenon view of accelerated system. 


A. Fermionic field quantization 

The Dirac equation in curved space-time is written by 


where i and j run over 1, 2, 3. By now, for squared 
equations we used / wcr , Xi and X 2 - 

To simplify these equations, we introduce the functions 
(j>± = Xi T X 2 and we can define and through 




^(x,w) 

^(x.w) 


(41) 


Given this definition, we can separate the equation for £ 
and £ in the form: 


(* 7 flV„ - m)^(x) = 0 , 


(32) 


where 


/ 1 2 \ 

X = (V, X ,X , U), 

w = (cj, k 1 , k 2 ), x=(x 1 ,x 2 ,u), 

(eo r = u~X, ( ei r = 8?, 
in = 

+ (33) 

In the Rindler coordinates, the equation becomes 

= (^y° mu ~ - iua l d^j (34) 

where a 1 = 7 ° 7 *. 

The fermionic field can be expanded as 


pOO pOO 


'k(x) = ^ / duj I d 2 k 
CT __j_ J 0 J —00 


h ,/,(+“) 4 . 

u wa Ywcr 1 rw—o 


We will the solution in the form 

^ = / W ff(x)e _i "^° 

for —00 < uj < 00 . 

The function / wcr is eigenstate of Hamiltonian as 


where 


H = a 


H /wcr — 


0 ■ ia 

mux - n - lUa Vi 


We denote two-component spinors Xa by 

f i'v'i = r xi( x j w) 
/ WCT (xj_ _ 2(X;W) 


(35) 

(36) 

(37) 

(38) 

(39) 


5 lJ udiiudj^) = m 2 u 2 + | iu ± - 


8 t 3 udi(udj£ ± ) = m 2 u 2 + ( iu =F 77 


^,(42a) 


(42b) 


Here, we introduce i as l 2 = (fc 1 ) 2 + ( k 2 ) 2 + to 2 by which 
the equation can be written as 


- ud 3 (ud 3 ) + ( iu ± i 


\ d 2 + d 2 2} C ± = J_ 

VT 

= -l 2 ^. 

The solution can be written in the form 

^ = A ± ('w)X{x 1 , x 2 , k 1 , k 2 )M±(u). 
This leads to the following decoupled equations 
[dl + d\- m 2 ]X(x 1 , x 2 , k 1 , k 2 ) = 0 , 


(43) 


(44) 


(45) 


d 2 1 d ( iw =F 7 ' 

+ — - 1 - 1 


d(£u) 2 lu d(£u) 


lu 


M±(u)=0. 


We can readily solve the equation in the form 

X(x 1 , x 2 , k 1 , k 2 ) = e zkaX , 
Mt{u) = K iQzf i(£u), 


(46) 


(47) 

(48) 


where the index a run over 1 and 2 . 

Using the arbitrary function A± and B ± of k, we find 
the most general solution can be written as 


£± = A ± (w)e ikaX K i , 


'■iu.(^ u )l 

<;±=B ± (w)e ik “* a K iQ±i (tu), 


(49a) 

(49b) 
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which implies 


The creation and annihilation operators should obey 


,/w (7 - 


1 

2 


c + +r i 
c + + r 
^ +r 
-c + + C" J 


r A+K^i^ + A-K^iu) 
,ik a x a B+K^i^u) + B~I\ iCj _i{iu) 
2 -A+K^itu) + A~I\ iCj+ i (tu) 
_-B+K^ +h {iu) + B~K iGl _ h {lu) 


.(50) 



It is easy to see that / W(T satisfies 


ia l di — mj° + 

\ 

ia-’dj — m'Y + 


( ic? \ 

1 


— 

V 2 ) 

U 


(io? \ 

l 


— 

V 2 ) 

U_ 


Ua = 0.(51) 


Therefore, the solutions of Dirac equation are simply 
given by 


J W O’ 


ia l di — mj° + 



(52) 


If we set 


A + = B- = 0 (53) 

we can find a normalized A~ and B + by letting 

1 

2 = N, (54) 

where 


A~ = B + = 


coshc 


7 t£ 


1 = 


e 

5 

a 


(55) 


which is normalized with respect to the inner product 



= S 3 (w - w')5 aa >, (56) 

where ij> = 7 0 and S is set to be v = const.. In this 

way, we find concrete form of fermionic field is written 

by 

p -iu)v/a-\-ik a x a 

^(x) = - —3 - u<?)(u, w), (57) 

(27r) 2 

where 


\u, w) 


= IV 7 0 


(ha 7° + rh)K iCj+ i(lu) + il'fK^Aiu) 


and 



- 1 - 


- 0 ' 


0 


1 

u+ = 

1 

, U- = 

0 


. 0 . 


. -1 . 


U(j 

(58) 


(59) 


B. Calculation of Cross Section 

The process of (3 decay in the accelerated frame looks 
very different from that in rest frame. In this case, a 
proton is stable but the whole space is FDU thermal bath 
characterized by a temperature T = a/2i:. Therefore, 
proton absorbs e~ and v from FDU thermal bath and do 
emit e + and v. Three processes are possible through the 
precesses: 


(i) p + + e“ —> n + u, 

(ii) p + + V —> n + e + , 

(iii) p + + e _ + V —> n. 


The transition rate is a combination of them. 

Formally, we can calculate the cross sections by same 
way of the inertial system but we have to deal with three 
processes. 

The transition amplitudes are written by 

= ( n \ ® |p) 

/*oo 

= — / dve iAmv (v ku(J ) 

a J -00 e e 

= ^S{u> e - -w v - Am)u^ ) 7 0 w^i" ) , (61a) 

f'TT v e 


~ Ue+ - A m)u^ e e + + \ 0 u^/\ (61b) 

«4(iiT) n = “ Am ) ft -”a7 ) 7°w^i” ) J (61c) 


where Si is replaced 7 ^ in inertial frame with 7 )) and 
u» = (a, 0 , 0 , 0 ). 

We assume that the observer is in the thermal bath. 
We attach the fermionic thermal factor for each process. 
Then the differential transition rate per absorbed and 
emitted particle energies for each processes are written 
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as 


1 


d^y 


T dui e - du>vd 2 k e - d?k v 


= l^r" n ^( w e-)[l-«F(w„)] 


CT _ (T v 


G% 


-(»») o K-) 

U<t„ TU a _ 


5{uj e - — (jj v — A to) 


2 5 7T 3 gTrAm cosh W<,-7T COSh d) v TT 


-,(62a) 


1 


d 6 ^ 


T cL> e -duJed 2 k e -d 2 kv 


= ^ mi Km" n F {u; e -)n F (u)v) 


T 


G 2 f 


-(—uip) n (“.-) 


8(u) e - + Us — Am) 


2 5 7r 3 e 7r ^ m cosliu)p-7rcosha;p7r 


■,(62c) 


where 


d 6 P ( P- 


T du>vdu) e +d 2 kvd 2 k e + 

= n F (uj p )[l-n F (uj e+ )\ 


CT e+ 


G% 


-( w e+) 0 (“») 

U,j l TUa p 


S(ujjy — u e + — Am) 


2 5 f 3 girAm cosh Ap 7 r cosh d; e + 7 T 


«f(w) = 


(63) 


is the fermionic thermal factor and T = 27rd(0) is total 
proper time of the proton. 

, (62b) We can obtain the completeness relations by 


Y w )u^\u, w) = A 2 7 ° 2? K iQ+ 1 (£) + *mi{(y 3 - 7°)A; 2 - + i (£) + (j 3 + J°)F 


-2 

x i£>— i 


(64) 


Using these relations, a direct calculation yields the spin sum for the process (i) as 


EEK 


(«„)-.o„K- ) 


7 


<T e - <X„ 

2 2 

= — r cosh u) e - 7r cosh a>^7r 

7T 4 


^ e _+i(4-)^„ + |(4)f +m e m,Re{A 2 , e +| (4-)A 2 -_ | (4)}] • (65) 


We can perform the analogous calculation for the processes (ii) and (iii) . 

According to Eq. (65), the differential transition rate of the process (i) is given by 


1 




T dw e -dw„d 2 fc e -d 2 4 


G 


-8(u> e - — u) u — Am) 


4-4 


K id e -+i (4- ) K iu, u +^ (4)| 2 + rh e rh „Re {^ e _ +i (4-)^„_i(4)} 


.( 66 ) 


By using them, the cross sections for each processes can be obtained as 


nr =f I Kr 

2 G 2 f 


^-7g7rAm 


/•OO 

r oo _ 

/. d<D e - 

/ d 2 fc e -4- 

J Am 

do 




*"«-+■ 


r (4- ) 


d 44 


7\" - , 

i{Q e - —A m) + 


i(4) 


+m e m 1/ Re<{ / d 2 k e -K 2 - _ + i(l e -) d 2 4A" 2 44) 


(67a) 


0/^2 /*c 

F I 

(") jTgiAni 7 0 


d<4 


d 2 4+4 


-^4zi + +i(4+) 


d 2 4 4 


K. 


i(oi + -Am)+5 ' 


+m e m„Re ^ j ^ d 4+dG-^ + + i (4+) J ^ d 4d4^ + _£^_i (4, 


, (67b) 
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pp—m z F 

(***) ^-7 g7rAm 


p Am 

pOO _ 

/ dw e - 

/ d 2 k e -l e - 

'o 

JO 


K iQ + i(i e ~) 


d 2 kali; 


K 


i{Co e - — Am)+-|' 


+m e ?n l/ Re< / d^kg-K 2 ^ ,i(£ e ~) / d 2 kpK 2 


i(uj e - —Am)— 2 


i(4) 


(67c) 


By summing them we can lump together the integration, and the total cross section becomes simply in the form 


pp—m _ -pp^n . pp >n . pp^i 
1 acc — 1 (i) ' 1 (ii) "T 1 (iii) 


7r' e" 


pOO 

/ d 2 k e i e 


2 r°° _ 

/ d 2 kj v 

Jo 


Jo 


A' 


i(tD-Am) + 


i(C) 


+m e m y Re^ / d 2 k e K 2 i(I e ) / d 2 k v I\ 


i(w-Am)- 


i(4) 


.( 68 ) 


It is hard to deal with these integral of modified Bessel’s. Of course, we can use Eq. (|2S|) again like the case of Section 
I. But we use more useful formula (see, for example, p219 of Ref. |e| or Ref. |fi~|). 


x°K v {x)K^{x) = ^Gf 4 x 2 


2 CT, 2 a + 2 


l(is + p, + cr), i(^ — fi + a), h(—v + /i + ct), |(—y — fj, + a) 


(69) 


and by using the definition of G function, the integrand can be represented by power of t and we can easily integrate 
with respect to t. 

We find 


r 2 

p p—>n 'Grp 


ds f dt 


dur 


2 3 7r 4 e TrAm J Cs 2-rri J Ct 2-rri J (2s + l)(2f + l)r(— s + l)r(— t + 1) 
r(—s + -)r(— t + -) r(—s + *(<d — Am) + i)r (—t + iu> + 1) 

+m e m l/ T(—s — — )r (—t — -)Re |r(— s + i(u> — Am) + l)r(— s — i{u> — Am))r(-f + iur + 1)E (—t — tw)| 


,(70) 


where all poles of complex s and t planes are picked up with C s and Ct, respectively, by definition of G function. 
To integrate with respect to w, we use the formula of Barnes [TL4|: 


/ ioo 

dur T(a + u)T(b + (D)r(c — A)r(d — ur) = 2-rri 

-ioo 


T(a + c)T(a + d)F(b + c)r(b + d) 
T(tt T b T c T d) 


[Re a, Re b, Re c,Re d > 0]. (71) 


Eventually, we find that the total cross section in the Rindler frame is 
" 5 Gp f ds f dt (ml) s (rh 2 y 


■pp —>n _ 

1 acc 


2 5 


7T2 e" 


>c s 


2-rri 


2-rri T(-s - t + 3)E(-s - t + |) 


r(-s -1 + iAm + 3) r(-s)r(-t)r(-s + 2)r(-t + 2 ) 


+Re |r(— s — t + iAm + 2)T(— s — t — iAm + 4) j T(— s + -)T (—t + -)T (—s + ^)T (—t + ^) 


■ (72) 


Comparing this to the results in inertial frame & we find that resulting expression agrees perfectly. This result 
shows the existence of Unruh effect is inevitable. ,_ 


IV. DISCUSSIONS 


analytic expressions for both frame and found that they 


We have analyzed the f3 decay and the inverse (3 de¬ 
cay of the accelerated proton in both frames. We found 
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Im s 



FIG. 1: All residues of F(— s) and r(— s — n) are picked up. 


The last form is the definition of modified Bessel function 
K^(z) for non-integer ft. And you find the formula in case 
/r is integer n by setting fi = n after partial differentiation 
of this formula by /r. 

Next we integrate Eq. ([l0|) as FIG.2. 

Im t 



agree with each other. If you see the calculation of G. 
E. A. Matsas and D. A. T. Vanzella in two-dimensional 
model 0, |ll|, you can realize that on four-dinrensional 
model it became hopelessly complicated integral. So the 
main problem in this time is the complication of inte¬ 
gral. To solve it, we used Barnes type representation as 
you can see from Eq. (0) and Eq. (|09|). We can demon¬ 
strate these formulae by picking up the poles in complex 
plane of the integral valuable. The sum of these read¬ 
ily becomes the infinity series which defines the special 
function. By using them, we accomplished perfectly an¬ 
alytical proof. 

It is straightforward to apply our technique for two 
dimensional setup used in Ref. 0 0, we can easily 
prove that the decay rates is independent of the frame. 


APPENDIX A 


FIG. 2: The contour separates the poles of T(—t) from those 
of T(f — z). 



dt T(-t)r(t-z) t+x t 

2ni T(— z) 


If you select the contour (i) you find 


,, V (~l)"r(-z + n) (B\ 
h n!T(-*) \Aj 



(A2) 


In this appendix A we proof both key equations 
Eq. ( p8| ) and Eq. (p0|). 

Firstly we derive Eq. (0). Through the path Ci, all 
poles of T(—s) and T(—s — /j) are picked up (see FIG.l). 


1 

2 


J C J S r( -* )r( -* - (|) 


£ \ 2 s~\~ fi 




2 n±/^ 


n =0 


1 00 

^EE 

n —0 ± 


(- 1 )" 


r! T(n ±/r + 1) sin(—n =F m) 71 " 


/ z \ 2n ±M 

V 2/ 


7T -J m (z)+/- m (z) 

2 sin p.7r 

(Al) 


So this integration is the expansion form of (A + B) z in 
B < A. Similarly if you select the contour (ii) then you 
obtain the expansion form of {A + B ) 2 in B > A. This 
integral representation of expansion of (A + B) z includes 
both cases of B < A and B > A by selecting the contour 
(i) and (ii), respectively. 


APPENDIX B 


In this appendix B, we show the explicit form of the 
integral. There are poles of the power of 1, 2 and 3. 

The integral can be simply calculated by the change of 
variables to 


s — i, t 


(Bl) 


where I^{z) is modified Bessel function of the first kind. After this transformation we obtain 
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pp 


—>n 


ci 5 Gp I" dt I" ds {drily t {rhl) t 

25 7r | e irAm J 2 iri Jq s 2 iri T(— s + 3 )r(— s + |) 


x 


r(-s + iAm + 3) 


2 T(-t)T(-t + 2)T(-s + t)T{-s + t + 2) 


+Re 


|r(— s + iAm + 2)r(— s - iAm + 4) j r(— t + -)T(— t + 1)T(— s + t + -)T(— s + 1 + 1) 


(B2) 


where if>{z) = In F(z) and C't is the path separating the poles of, for example in first term, T(— t)T{—t + 2) from 
those of r(-s + t)F(—s + t + 2), and C s is the path which picks up all the poles in s complex plane. 

Firstly, we start by t integration because the existence of poles in t complex plane is independent of s. If we 
integrate by t, we obtain for m e > rh v : 


r*»" = ■»*<%_ f *- Al - J ( r( -»)r(-» + 2 ) + (A ) r(-» + i)r(-« + 3) 

2 5 7T2 i 27tj F(—s + 3)r(— s + ^ ) 1 l \ m e 


OO / ~ 

+ V ( ^ 


2 n r 


tp{n — 1) + y(n + 1) — ijj{—s + n) — ip{—s + n + 2) — 2 In 




r(— s + n)r(-s + n + 2) 
r(n- l)r(n+ 1 ) 


r(-s + iAm + 3) 


E t m^ 

' 77 l e 


n —0 


2n+3 


ip(n + 1) + ip(n + 2) — ip{—s + n + 2) — ij){—s + n + 3) — 2 In 


m e 


xRe |r(— s + iAm + 2)r(— s — iAm. + 4) j 


r(— s + n + 2)r(— s + n + 3) 
r(n + i)r(?z + 2 ) 


(B3) 


The simplest case is for massless neutrino. If we set m„ = 0, then only the first term remains non-vanishing. After 
a valuable transformation s —> s + | we have 


r p^n 
1 m„ =o 


a 5 G 2 F rhl 

I' ds 

T (—s + iAm + §) 

2 5 7 r 5 e 7rAm J 

2ni 

B(-s + | 


—mi 


}m e a‘ 

7 — 24 l o - 

327T2 e 7rAm V a 


3 

2 ’ 


2 a 


(B4) 


This is exactly the cross section obtained by Vanzella and Matsas with cy 


1 and ca = 0 (see Eq. (4.19) in [p~5|) . 


We perform the integration of Eq. (B3) with respect to s and find the following expansions of the cross section 
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pp —>n _ 


EE 

. m—0 =h 


D-h^' + E 


i/j(n — 1) + ip{n + 1) — ii>(n — m =F iAm — 1) — i/j(n — m =F * Am — 3) — 2 In 


n —0 n —2 

(—l) m r(n — to =F iAm — l)T(n — to =F iAm — 3) / 


T(n-l)T(n + l) 


x ■ 


(—m 2iAm)B(m + 1, —2 m =F 2iAm) 

4 


V 1 


— (1 — 7) r(«Am+l) 


m v 

~2 


m e 

2 


^ r(2fc - 5) r(-fc + iAm + 3) 

fc= 3 


r(fc-i)r(fc + i) 


W t 

T 


2k 


2k 


°° r(2fc-3) r(—fc + *Am + 2 ) 

fc=2 


r(fc- i)r(fc + i) 


TOi 

T 


2k / - \ 2 


+ 


m h 

~2 


2 / ~ \ 2 k 

TO, 


OO OO 


+EE 


—8 - i/’(2[n + to — 3] + 1) + 2 ^ — mi iAm + 3) + E V>(n ± 1) + ijj(m ± 1) — 2 In 


TO,,TO e 


r(2[n + to — 3] + 1) r (—n — m . + *Am + 3) 


r(n - l)T(n i l)r(w - l)r(m + 1 ) 


m L 

~2 


2 n / ~ \ 2m 

m e 


OO OO 


+EE-E 


n—0 m=0 


il>(n + 1) + + 2) — if>(n — m -F iAm) — i/j(n — to =p iAm + 1) — 2 In- 

m e 


E 


(—l) m r(n — m =F *Am)r(n — to =F iAm- + 1) / fn v \ 

2 2 T(rc + l)T(n + 2 )(—to =p 2iAm + 2 )B(m + 1, —2?n 2iAm + 2) V/ 

ijj(n + 1) + if>(n + 2) — ip(n — m -F iAm — 2) — ip{n — m -f iAm — 1) — 2 In -A 


2n+3 


2(m±iAm+l) 


m e 


(—l) m r(?z — m ± iAm — 2)r(n — mi iAm — 1) 


:) 


2n+3 


2(m=piAm+4) 


2 4 T(rc + l)T(n + 2)(—m ± 2iAm — 2)B(m + 1, —2m ± 2iAm — 2) \ TO e 
—8V’(2[n + m] + 1) + 2 i/j(—n — mi iAm q= 1) + ip(n i 1) + ip(n i 2) + ij)(m + 1) + ip(m i 2) — 2 In mv ™ e 


T(2[n i m] + l)r(— n — to + iAm — l)r(— n — m — iAm + 1) /TOj, 


T(n + l)T(ra + 2 )T(to + l)r(m + 2 ) 


r 


V 2 


2n+3 / ~ \ 2m+3'\ 

f m e 

V 2 


(B5) 


where B(p, q ) = and 7 is Euler’s constant. 

This is the final form of the cross section. A natu¬ 
ral question is that we have obtained the result which just by interchanging rh u and m e for m e < m v . This can 


metric. The resolution of this puzzle is that the integral is 
of discontinuous type. Namely, we can obtain the result 


is not manifestly symmetric with respect to rh u and m e 
although the original expression (pTl) is manifestly sym- 


be checked directly by changing the order of integrations. 
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